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On the quantization accuracy of the acoustoelectric current 
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By deriving analytical formulae for the quantization accuracy of the acoustoelectric current, we 
reveal that: 1) the flatness of the current plateau for the typical present devices has reached the 
theoretical limit of about lOOppm over a 1/1000 change of the gate voltage; 2) increasing the 
transport channel length, and counterintuitively, increasing the acoustic wavelength as well, would 
improve the quantization accuracy, and very promisingly up to O.Olppm as required for a quantum 
current standard. 
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Introduction: The observation of a quantized electric 
current driven by surface acoustic wave (SAW) in 1996 [1| 
visualizes a possible quantum current standard. The de- 
vice comprises a split gate build on a GaAs-AlGaAs het- 
erostructure containing a two-dimensional electron gas 
(2DEG). A negative voltage Vg applied to the split gate 
depletes the 2DEG beneath the gate and leaves a one- 
dimensional (ID) transport channel. When pinched off, 
the ID channel becomes a saddle-like potential barrier. 
At large wave amplitude, the SAW-induced potential 
minima become moving quantum dots (QDs), which cap- 
ture electrons within the dots and drag them through the 
potential barrier. The important finding of Ref. [l|] is 
that such acoustoelectric current displays plateaux over 
certain range of Vg and/or SAW power, and the plateaux 
fit well into the quantized expression / = NefsAW, with 
A'' an integer, e the electron charge, and /saw the SAW 
frequency in the GHz band. 

Qualitatively, this quantization phenomenon is ex- 
plained by the Coulomb-blockade effect: At a certain 
depth the QD can only hold a fixed number of elec- 
trons, because excess occupation is prohibited by signifi- 
cant Coulomb repulsion among electrons within the small 
QD. The detailed electron transport process and possi- 
ble error mechanism in such devices have been studied 
theoretically but quite controversially, in both classical 
dynamics [2] and quantum mechanism [S, U, a IgI, LZf. 
Sometimes, a very accurate current quantization (bet- 
ter than 10~^°) is predicted [g|. However, the flatness 
of the current plateau achieved in the o rig inal [l| and 
subsequent experiments [3, |9|, [lOJ, |ll|, [l^, |l3| never goes 
beyond ~ 100ppm/10~^V"g, far below the level of ~ 10~^ 
required in metrology. 

In this paper, we aim at a clear clue on what factors 
determine the quantization accuracy of the acoustoelec- 
tric current. We start with a novel analysis of the elec- 
tron transport process and possible error mechanisms by 
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clarifying on the different mechanisms for electrons to 
leave a QD. Then by some approximation techniques, we 
derive analytical formulae for the quantization accuracy. 
These formulae reveal clearly the limitation of the typical 
present devices and the promising possibility of improv- 
ing the accuracy for metrological use. 

The electron transport process and error mechanisms: 
Since the quantization of acoustoelectric current is essen- 
tially the quantization of electron number in each QD, we 
find it very illuminating to first study how the electrons 
can escape a QD. We clarify that there are three such 
escape mechanisms when the QD climbs up the potential 
barrier: 

(i) directly fiowing out when the QD reduces its depth; 

(ii) back tunneling through the rear wall of the QD; 

(iii) non-adiabatic transition into excited states, fol- 
lowed by an easier tunneling. 

Mechanism (i) occurs because the depth of the QD de- 
pends on the slope of the potential barrier, which is not 
uniform and typically has a maximum at the waist of 
the barrier. For mechanism (iii), the excited state can in 
principle be a free state so that the electron escapes the 
QD directly, but we focus on the dominant case of tran- 
siting to a higher-energy bound state in the QD, and the 
electron escapes via an enhanced tunneling probability. 

With the above clarification, we divide the electron 
transport process through the potential barrier into three 
segments, characterized by the location of the QD: 

(a) before the steepest (S) point of the barrier; 

(b) from the S point to the top of the barrier; 

(c) after the top of the barrier. 

In segment (a), all three escape mechanisms are in ef- 
fect. The QD becomes shallowest at the S point, so in 
segment (b) electrons can no longer fiow out directly and 
only mechanisms (ii) and (iii) remain. Segment (c) is ir- 
relevant for us because all escape mechanisms effectively 
kick the electrons into the drain instead of the source. 



Let us denote Nx as the QD occupation number at the 
entrance of the ID channel, Ns as that at the S point, 
and Nt as that at the top of the barrier. As we explained 
above, Nj- can be take as the final number of transported 
electrons per QD, so the acoustoelectric current is 



/ = NrefsA' 



w- 



(1) 



Variation in Nt means an error in the current quantiza- 
tion. It should be noted that Nt can only differ slightly 
from Ng , by the "mild" escape mechanisms (ii) and (iii) . 
In contrast, Ng may differ dramatically from Nx, be- 
cause in segment (a) the QD reduces its depth greatly, 
so that electrons can escape in a "wild" manner of flowing 
out directly. We write 



NT^il-Pt~Pv)Ns. 



(2) 



Here Pt is the (small) probability for an electron to escape 
the QD in segment (b) by direct tunneling, and Py is 
the (small) probability for an electron to first transit to 
excited states by non-adiabatic correction, then escape 
the QD by an easier tunneling due to higher energy. From 
Eq. ([2]), we can define two additive error mechanisms in 
the quantization of / = Ntc/saw- First, the "wild" 
reduction from Nx to Ng could be inaccurate, which 
results in an indefinite Ng- Second, even if Ng = 1, viz., 
exactly one electron remains in each QD at the S point, 
this single electron may still escape the QD in segment 
(b) by non-zero Pt and Py. 

Investigation of the "wild" reduction from Nx to Ng 
necessarily involves many-electron dynamics. It has been 
calculated numerically with classical dynamics in Ref. 
[21 . And quantum mechanically, Ref. f5] has performed 
a phenomenological study of this process at the very en- 
trance to the ID channel. In the following, we concen- 
trate on the "mild" reduction from Ng to Nt, starting 
with Ng = 1, which is then just a one-particle problem. 
Experimentally, the first current plateau (corresponding 
to one electron transfer per SAW cycle) usually has the 
best quality. As we will see below, the escape probabil- 
ity {Pt -t- Py) alone explains the observed quantization 
accuracy very well. This means that the error in the re- 
duction from Nx to Ng is at most of the same order as 
(Pt + Pv). 

Back tunneling of electrons out of a moving QD has 
been studied numerically in Refs. [3|, y, |7|. Here, we de- 
velop some approximation methods which allow us to es- 
timate this effect analytically. Let us first look at the tun- 
neling probability given by the familiar WKB formula: 



P 



exp 



dxy/2m*{V{x)-E) 



(3) 



Here m* — O.OGTrrie is the effective electron mass, V{x) 
describes the QD rear wall to penetrate, E is the electron 
energy, and a;i.2 are the two transition points at which 



V{x) — E. It should be noted that P in Eq. ^ is derived 
as a transmission coefficient, viz., the intensity ratio of 
the out-going and incident beams. In other words, P is 
the penetration probability per collision. If an electron 
captured in a QD collides with the QD rear wall F times 
per unit time, then the total probability for it to tunnel 
out of the QD in a time duration [ti , t2\ is 



Pt = 



PTdt. 



(4) 



The collision frequency F is just the electron oscillation 
frequency in the QD. If we approximate the electron mo- 
tion in the QD as a harmonic oscillator with angular 
frequency Wg, then T — ff, = tt'e/27r, which is typically 
much higher than the SAW frequency /saw- Evidently, 
Eq. (JD) is equally valid if P and F vary with time. 

Eq. ^ tells us that the back-tunneling probability P 
is largely determined by the height and width of effec- 
tive part of the QD rear wall lying above the electron 
energy. To examine these parameters, we make a reason- 
able approximation by considering the superposition of a 
sinusoidal potential wave Asin(27ra;/A) and a straight po- 
tential slope Kx, with k corresponding to the local slope 
of the real potential barrier. The largest derivative of 
the sinusoidal wave, 2'kA/X, must exceeds k, otherwise 
no QD can be formed. In a recent paper [ij], the SAW- 
induced wave amplitude A was estimated to be 25 meV. 
In contrast, to make iV^ = 1 the depth A of the QD at 
the S point must be smaller than the chargin g en ergy of 
the QD, which is estimated to be ~ ImeV [a llJI- Since 
now A <C A, the slope kq at the S point must be very 
close to 2itA/\ so as to produce a superposed QD much 
shallower than A. We therefore define 

with 77 a small parameter. Then, a slight algebra shows 
that the depth A of the QD is approximated by 

3 



^/2 



i^^A, 



(6) 



and the effective width of the QD rear wall is approxi 
mated by 

3V2 



W ■ 



2-K 



■^\. 



(7) 



By symmetry, W is also the effective diameter of the QD. 
We see that at the S point of the potential barrier, the 
size of the QD is much smaller than the SAW wavelength. 
The parameter 77 is not free. It can be fixed as follows: 
To permit only one electron in the QD, its depth A must 
not exceed the charging energy E^ ^ e^ /rteS, with r = 
W/2 the effective radius of the QD and ecff the effective 
dielectric constant in the QD. From A < £■£, we have 



< 



Ait e^ 



9AA 



toff 



(8) 



The upper limit corresponds to the best quantization 
accuracy. To have a quantitative feeling of these pa- 
rameters, take A — 25meV, A = 1/im, and tentatively 
Eoff — 13, we find 

7] < 0.079, A < 1.2meV, W < O-ig^im. (9) 

They are all consistent with our preassumptions. 

There is a hidden subtlety in Eq. ([8]) which calls for 
special attention: the choice of e^s- Quite often, it is 
argued that since the 2DEG resides very close to the 
surface, (viz., with a distance d ^ 0.1/im much smaller 
than the SAW wavelength) eeff is roughly half of the value 
for a body material. However, Eq. ([9]) indicates that 
the size of the QD is also much smaller than the SAW 
wavelength, and actually comparable to d, thus within 
the QD Ccff can be close to the value (~ 13) for a body 
material. 

To apply Eq. ((3]), we still need the electron energy 
inside the QD. We approximate the electron motion as a 
harmonic oscillator with an angular frequency LOe' 



* 2 



Kq 



■^(¥>V'-<2,^/A' 
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Ai^f^. (10) 



A quick ordcr-of-magnitudc estimation of the tunneling 
probability can be obtained by approximating the QD 
rear wall as a triangle with height A and width W: 



exp 



2 2 A- 



E 



W)^/2m*{A-E) 



(11) 



Here E ~ {n + ^)fiLOf. is the electron energy relative to 
the bottom of the QD, and the ground state has n — Q. 
The final ingredient needed to compute Pt is an explicit 
expression of the potential barrier. We take it to be 



V{x) ^Va/cosh^ix/L). 



fl2) 



The effective width of this potential barrier is 2L, the 
steepest point is at xq ~ — 0.66L, with a slope 



Ko = Ciyo/i-0.77Vo/L. 



(13) 



This is our first master equation. Here Ct and Ce are 
straightforwardly computed coefficients independent of 
A, L, and Vq. We postpone discussing the impfications 
of Eq. (jlSp . and turn to the other error mechanism: 

Non-adiabatic excitation: This effect can be most con- 
veniently studied by adopting a coordinate system co- 
moving with the SAW. Because the potential barrier is 
not uniform, even in this co-moving coordinate system 
the electron feels a varying external field, which can lead 
to quantum transitions. Namely, an electron in state a 
at time t = has a non-zero probability Pba(t) to be in 
another state 6 at a later time t. To first order, Pba{t) is 
given by [15| 



Pba{t) 



^^-^^ba' 



dH \ 



sm^[uJbat/2] 



(16) 



Here ujba{t) — [Eb{t) — Ea{t)]/h, with Ea.b{t) eigenvalues 
of the instantaneous Hamiltonian H(t). (^)i is the 
matrix element of ^ between the states b and a. Be- 
cause the electron oscillation frequency is much higher 
than the SAW frequency, the factor s\T?[uobat/%, in Eq. 
(I16p oscillates very quickly. By approximating the elec- 
tron motion in the QD as a harmonic oscillator, our prob- 
lem resembles the textbook example of a charged har- 
monic oscillator in a time-dependent electric field [l5|, 
for which only the transition into the first excited state 
is relevant, and the transition probability Pio can be ex- 
pressed as 



P^ 



10 



1 d^V 



m*uj'? 



dx^ 



sm'^[uJet/2]. 



(17) 



Here v,, = ^2?iuje/'m* is an effective electron speed in 
the QD, with uje from Eq. PH)) . To see the dependence 
of Pio on the experimental parameters, we write 



d^V/dx^ = C-iVo/L^ 



(18) 



with the coefficient C2 independent of Vq and L. Follow- 
ing the same steps as in deriving Eq. (|15p . Pio can be 
expressed as 



The total probability for an electron to tunnel out in 
segment (b) is then 



Pt = 



P{x) 



UJe dx 

2tt v^ ' 



(14) 



where Vg — 2800m/s is the SAW speed. 

There are three free parameters in the above expres- 
sions: the SAW wavelength A, the barrier width 2L, and 
the barrier height Vg. From Eqs. ([ MT5)) . the tunneling 
probability at the S point can be expressed as 



P ~ exp 



-CtXiL-^V^ 



l-CA--L-^V'^ 



(15) 



Pio ~ C,X-L~-Vo^ sm^[uj,t/2] 



(19) 



This is our second master equation. Here Cy is another 
coefficient independent of A, i, and Vq- 

Quantization accuracy for the present devices: We now 
plug in explicit parameters to estimate the theoretical 
limit of the quantization accuracyof the acoustoelectric 
current. In a typical experiment [l|, [sj, |9|, [lOJ, [llj, ll2l.ll3J. 
one has A = 2i = l/ini. There is just one parameter 
to tune: the gate voltage (which determines the barrier 
height Vq) or the SAW amplitude. But this parameter 
does not bring much freedom either. Vq must be at least 
several tens of meV in order to be well above the electron 
Fermi energy. Furthermore, our two master equations 



(115)) and (|19p show a rather week dependence on Vq- For 
the experimentaUy estimated A = 25meV IJ], Eqs. ^ 
and (fT^ give Vq ~ 94nieV, in consistent with our ex- 
pectation. Then, Eq. pT|) or ([TS]) gives the tunnehng 
probabihty of the ground-state electron at the S point: 
Pq ~ 10^^, and Eq. pi]) gives the total probability for 
this electron to tunnel out in segment (b): Pt ~ I.SPq. 
This integrated tunneling probability is of the same or- 
der as the tunneling probability at the S point because 
P drops very quickly as x moves away from the S point. 

The situation for non-adiabatic transition requires 
some caution. For an electron to eventually escape the 
QD by this mechanism, the transition must be followed 
by an easy tunneling. The tunneling probability is largest 
at the S point, but Eq. (fTT]) says that the transition prob- 
ability is related to the second derivative dV'^ /dx^, which 
is zero at the S point. At first thought this seems to imply 
that the non-adiabatic mechanism is not relevant. How- 
ever, we note that our whole discussion for non-adiabatic 
transition can as well be applied to the center-of-mass 
motion of the interacting electrons in segment (a) , where 
d^V/dx^ can be significant. This means that during this 
period the center-of-mass motion of electrons may tran- 
sit to an excited state, by a probability of the same order 
as given by Eq. (fT7|) . The excited electrons will then es- 
cape the QD by the significant tunneling probability near 
the S point, leading to an error in current quantization. 
Therefore, to safely control this kind of quantization er- 
ror, we should estimate P^ as the product of the maxi- 
mum Pio in segment (a) with Pi at the S point, where 
Pi is the tunneling probability of the electron in the first 
excited state. A slight calculation with Eq. (fT7)) gives 
Pio -- 2 X 10~*, and Eq. HID gives Pi - 0.5. Note that 
the WKB formulae is only valid at very small tunneling 
probability, therefore a Pi ~ 0.5 should be effectively 
taken as one (viz., 100% electrons escape). So we finally 
estimate P„ in Eq. ([21) to be ~ 2 x lO""*. 

Putting Pt and P^ together, we see that the acous- 
toelectric current at the first plateau deviates from the 
quantized value /i = e/sAW by at least 1.5 x 10"'^. This 
absolute deviation is hard to detect due to the zero-point 
drift in measuring the electric current. What is easier 
to detect is the slope of the current plateau. By fixing 
A = 25meV and varying Vq, the slope is computed to 
be dl/dVo ^ lO^'^/i/lO^^Vb- Since Vq is proportional to 
the gate voltage Vg, this means that the acoustoelectric 
current varies by lOOppm as Vg varies by 1/1000. This 
agrees very well with the best experimental result [9| . 

Tuning the parameters for better accuracy: From our 
two master equations (|15l) and ([19]), we see that the quan- 
tization accuracy of acoustoelectric current is essentially 
set by the device parameters A and L. A longer transport 
channel, viz, a wider potential barrier, always reduces 
the error in current quantization. And quite counterin- 
tuitively, a longer SAW wavelength may also be advan- 



tageous, because as A increases the tunneling probabil- 
ity falls exponentially while the non-adiabatic transition 
probability increases in a power law. 

A larger A also brings a hidden benefit that the QD 
size may indeed considerably exceed d, the distance of the 
2DEG to the surface, thus the effective dielectric constant 
reduces and the Coulomb-blockade effect ehances. E.g., 
if we take X = L = 5/im, Vq = 90meV, and pre-assume 
that r = iW^ <C d so that Ccff — 13/2, we get r ~ 0.46^m, 
in consistent with the preassumption. Proceeding with 
Eoff ^ 13/2, Eqs. mi) and dH]) give P^ ^ Pt ^ IQ-^^ 
and Pi - IQ-i", Eq. dTT]) gives Pio - 10-^ and thus 
Py ~ lO^^'^. The slope of the current plateau is computed 
to be dl/dV ~ 10~^^/i/10~^Vb. These accuracies are far 
enough for metrological use, which requires O.Olppm. 

Summary: We performed a novel examination of the 
quantization accuracy of the acoustoelectric current. By 
clarifying the electron transport process and error mech- 
anisms, we obtain a clear clue on the theoretical limit of 
the quantization accuracy, expressed by several analyt- 
ical approximate formulae. We have shown unambigu- 
ously that the so far achieved quantization accuracy is 
already the theoretical limit of the present devices with 
A = 2L ~ 1/zm. To improve the accuracy to meet metro- 
logical requirement, it is necessary, and very promising, 
to have a larger L and/or A, e.g., of 5/xm. 
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